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Abstract 

Corner transfer matrices are a useful tool in the statistical mechan- 
ics of simple two-dimensinal models. They can be very effective way of 
obtaining series expansions of unsolved models, and of calculating the 
order parameters of solved ones. Here we review these features and 
discuss the reason why the method fails to give the order parameter 
of the chiral Potts model. 
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1 Introduction 

In an early work of the author, [1] he evaluated numerically a sequence of 
approximations for the free energy (or entropy) of the monomer-dimer system 
on the square lattice, using a variational approximation for the eigenvalues 
and eigenvectors of the row-to-row transfer matrix. His interest was to see if 
the model exhibited a phase transition as the density of dimers increased from 
zero to the close packed value of one-half. No such transition was observed: 
the model is critical only at close packing, where it can be solved exactly[20j. 
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This was reflected in the poorer convergence of the approximations as close 
packing is approached. 

The result may have been negative, but the work did lead naturally to 
the development of the useful concept of "corner transfer matrices" . 

These matrices build up the lattice by rotations about the centre. For 
an infinite lattice they are of infinite dimension, but for the zero-field Ising 
model, one can write them in terms of the spinor operators (or Clifford 
algebra) introduced by Kaufman [2T] and diagonalize them exactly. [21 [25] The 
results were surprisingly simple, and provided a reasonably direct way of 
obtaining the spontaneous magnetization (the order parameter) of the Ising 
model, compared with the previous verifications of the Onsager-Kaufman 
result [23] by Yang,[2S] Montroll, Potts and Ward[22] and others. 

It was also realised that the corner transfer matrices satisfy certain equa- 
tions. Again, to give exact results the matrices and equations must be 
infinite-dimensional, but a finite truncation corresponds to a variational ap- 
proximation of the type used for the monomer-dimer system. Such trun- 
cations can give useful approximations to the exact results, even for models 
which have not been solved exactly. For instance, Baxter and Enting[12] were 
able to obtain 24 terms in the low-temperature series expansion of the Ising 
model in a field, using only 15 by 15 matrices. This was a great improvement 
on the 12-term series expansions previously obtained. |24j 

Here we present these equations for a quite general "edge-interaction" 
model, where spins live on sites of the square lattice and the interactions 
are between adjacent spins. We discuss the simplifications that arise for a 
"solvable" model, i.e. one whose Boltzmann weights satisfy the star-triangle 
relation. We emphasize the importance of the "rapidity-difference" property 
and indicate how it ensures that the corner transfer matrices commute and 
have a very simple eigenvalue spectrum. This makes it easy to obtain the 
order parameters of such models. 

Finally we discuss the chiral Potts model, and indicate how the lack of a 
rapidity difference property for that model prevents its solution by the same 
corner transfer matrix techniques. 

For definiteness, here we consider only Z^- invariant models (such as the 
Ising, self-dual Potts and chiral Potts models), where each spin can take N 
values and the Boltzmann weights depend only on the difference (modulo N) 
of the spin values. The extension to more general edge-interaction models 
(perhaps including site weight functions as in the solvable Kashiwara-Miwa 
modeipiJ] - [T5]) is straightforward. 
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2 Square lattice edge-interaction models 



In Figure [T] we have drawn the square lattice C diagonally, denoting sites by 
circles and edges by solid lines. On each site i there is a "spin" <jj, taking 
some discrete set of N values, say 0, 1, . . . , N — 1. Spins on adjacent sites i 
and j, with j above and to the right of i, interact with Boltzmann weight 
function W(<7i — Oj), as indicated. Similarly, spins on sites k and I, with I 
above and to the left of k, interact with weight function W(ak — cri). The 
partition function is 

2 = En^-^on^K-^) > (2-i) 

the sum being over all values of all the spins on the lattice, the first product 
over all SW - NE edges, the second over all SE - NW edges. 
If the lattice has L sites, we expect the limit 

k = lim Z 1/L (2.2) 

to exist, and to be independent of the manner in which £ becomes large, so 
long as it becomes infinite in all directions. 

In Figure Q] we have singled out a central site with spin a. The average 
of any function f(a) is defined as 

(f(a)) = r^/^n^-^n^-^) • (2-3) 

a i,j k,l 

We also expect any such average to tend to a limit, provided we impose 
suitable conditions on the values of the boundary spins, and take a limit 
where a becomes infinitely far from any boundary. 



2.1 Corner transfer matrices 

Here we define the corner transfer matrices (and the associated row and 
column matrices) for a general anisotropic edge-interaction model. 

We have again drawn the diagonal square lattice in Fig. [2] and have 
singled out a central edge linking two sites with spins a and b. We have 
then used dotted lines to divide the lattice into nine parts: the central edge, 
four "corners" labelled respectively A 1 , B2, S4, and four columns or rows 
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Fi gure 1: The square lattice (solid lines, drawn diagonally), showing a central spin a 
and the (dotted) rapidity lines. 

Take all the spins on the outermost boundary of the lattice to have some 
fixed valued say zero. Let A be the set of spins on the dotted line to the left 
of the spin a and \x the set of spins on the dotted line below a. Consider all 
the edges of the lattice below A and to the left of //, i.e. all the edges in the 
lower-left corner of the lattice. Take the product over all such edges of the 
appropriate weight functions W or W, and sum over all the spins internal to 
that corner. The resulting partition function of the corner will be a function 
of a, A and \i. Write it as (A 1 ) AjA1 . 

Similarly, let (F^)^ v be the product of the weight functions of the edges 
between /i and v. Let (B 2 ) u ,t be the partition function of the lower-right 
corner of the lattice, i.e. the portion to the right of /i and below r, and so 
on. Then the partition function of the complete lattice is Z = Zyy, where 

Z w = J2W(a-b)(A*)^(F?)^ v (B 2 ) VjT ...(Gl) p , x , (2.4) 

^or an ordered system, the boundary spins should be set to values consistent with the 
ground state of the particular phase under consideration. 
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Figure 2: The parts of the lattice corresponding to the various matrices A,B,F, G. 

the sum being over all values of a, b and all the spins on the dotted lines of 
Fig. El 

Plainly we can regard (Aj)^ as the element (A, /i) of a matrix A X) (F 1 ) fJLyV 
as the element (fi, v) of a matrix F l , etc. Then (12. 4p simplifies to 

Z w = Y,W{a -b) Trace A a 1 F^B 2 G b 2 A b 3 F b 3 B A G a 4 , (2.5) 

a.b 

the sum now being only over the spins a, b. 

There are other ways of building up Z using the matrices Ai,Bi,Fi, Gi as 
building blocks, most of them being simpler than (I2.5p . They are 

Z w = J2 W ( a ~ h ) T™ 06 B^A^F* B z g\a\fI , 

a.b 



Z A = Tr ace A^A^Al , Z B = Trace B^B^ , 

a 



Zi = Trace A> i F i B i+i B i+2 G i+2 A i+3 , 

a 

where i — 1, ... ,4 and the suffixes containing i are to be interpreted modulo 
4. We have indicated how Za and Z 2 are constructed in Fig. [31 

Note that F i always follows A i and precedes B i+1 , while G i follows Bi 
and precedes A a i+l . 
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Figure 3: The partition functions Za, Z 2 . We have omitted the edges of £, 
but their positions can be deduced from the full circle denoting the lattice 
spin a. 

If we keep the sizes of the blocks corresponding to the A iy Bi, Fi, Gi matri- 
ces fixed (and mutually consistent), then Za, Z b correspond to the smallest 
lattices, Zw, Z^ to the largest, and Z\,...Z± to intermediate ones. Using 
( 12.21) . we find that in the limit of the lattices becoming large, the partition 
function per site is 

k = ZaZbZwZ W . (2 . 6) 
Z\Z 2 Zj,Z^ 

We shall find it convenient to define four Boltzmann weight functions 
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U 1 ,...,U i by 



U^n) = W{n) , U 2 {n) = W{n) , 

U 3 (n) = W(-n) , U 4 (n)=W(-n) (2.7) 

and to set 

Z\ = Z 3 = Z w , Z 2 = Z 4 = Zyy . (2.8) 

The equations [3], H] for the corner transfer matrices can be obtained by 
requiring that the expression (12. 6p be stationary with repect to variations 
in the matrices Ai, Bi, F^, G{. This variational principle leads to the four 
equations 



J2Ui(a- b)G i+1 A i+2 F- +2 B i+3 G a i+3 = rj { B i+2 G a i+2 A a i+3 , (2.9a) 



E Ui-i(a - tyFi^GUliF^ = A:f:b 1+1 , (2.9b) 



Fi B i+1 B l+2 G i+2 = £i A i+l A i+2 , (2.9c) 

E GUA-iAFt = & Bi-xB, . (2.9d) 
b 

Our notation is now such that rotating the lattice through 90° is equiva- 
lent to incrementing i by 1 in all suffixes. 

Pre-multiplying both sides of (I2.9al) by A^F^Bi+i and summing over a, 
we find that 

Vt = Z,/Z t . (2.10) 

Similarly, 

rt'i = Z l+1 /Z t , i i = Z l jZ A , & = Zi/Z B (2.11) 
and, for instance, 

« = ^V(& +1 £ +3 ) • (2-12) 

The equations (I2.9al) - (12.9dj) can be represented graphically, each side 
being a semi-lattice, e.g. for i = 3 the equation (I2.9al) can be represented as 
in Fig. H 

We can only expect these equations to be exact (for an infinite system) 
when the matrices are infinitely dimensional. However, if we take them to 
finite-dimensional, then the fact that fl2.9al) - f!2.9dl) are derived from the 
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Figure 4: Graphical representation of the elements (p, p') of the two sides of equation 
(j^a|) for i = 3. 

variational principle for n ensures that these equations are mutually consis- 
tent and that they will define the Ai,Bi,F{, Gi (to within irrelevant similarity 
transformations). There will be many solutions - a test of the utility of any 
particular solution is that it should reproduce as many of the largest eigen- 
values of of AiA^A^A^ and -B1-B2-B3-B4 as possible for that truncation. 

Such truncations can provide a powerful tool for performing numerical 
or series expansion calculations on otherwise unsolved models. [H] Indeed, it 
was such a calculation that led to the exact solution of the hard hexagon 
model. [HE] 

3 Solvable models 

We now focus our attention on models that are known to be solvable and 
satisfy the star-triangle relation. Their corner transfer matrix have spe- 
cial simplifying properties. These properties are true only in the infinite- 
dimensional limit, but for any given truncation they will be approximately 
true, in particular they will hold to appropriate orders in a low-temperature 
series expansion. 



8 



In general, if a model satisfies the star-triangle relation, then its weights 
W(n),W(n) depend on two variables p and q, called "rapidities". The vari- 
able p is assocaited with the vertical direction of the lattice C, q with the 
horizontal direction. We display this dependence by writing the weights as 
Wpq(n),Wpq(n). The star-triangle relation relates three sets of Boltzmann 
weights, with rapidities (p, q), (p,r), (q,r). It is 

£ W qr (b - d) W pr (a - d) W pq (d - c) 

d 

= R pqr W pq (a-b)W(b-c)W qr (a-c) , (3.13) 

for all p, q, r and values of the three external spins a, b, c. The sum is over all 
values of the internal spin d. 

A second relation must also hold, where all the spin-difference arguments 
in (I3.13P are negated, so that b — d becomes d — b, a — b becomes b — a, etc. 
Provided we impose cyclic boundary conditions, (13.131) ensures that the row- 
to- row transfer matrices of two models, one with rapidities {p,q), the other 
with (p,r), commute. [4] Thus it is natural to allow q, r to vary from row to 
row of the lattice C, and in Fig. [1] we have used this freedom to associate 
two rapidities q, q' with the horizontal rows of the lattice. Moreover, this is 
true even if the rapidities on the dotted vertical lines of Fig. [I] are different, 
e.g. if the lines to the left of a have rapidity p, while those to the right have 
rapidity p',as shown. 

More generally, we can allow a different rapidities p for each of the dotted 
vertical lines in Fig. [T] , and a different rapidity q for each horizontal line. 
The Boltzmann weight of any edge of C is W pq {n) or W pq (n), where the p 
and q are the rapidities of the dotted lines intersecting that edge. 

It follows that the eigenvectors of these transfer matrices are independent 
of horizontal rapidities such as the q, q' shown in Fig. [TJ The correlations 
of spins lying within a single row of £ will therefore also be independent of 
q,q'. 

But (to within an overall irrelevant normalization factor ) the expressions 
represented by Fig. H] are precisely such correlations, involving only the spin 
sets p, p' and the spin a on the upper edge of each side. This implies that 

(BiGiA^p^ = independent of q, q' . (3-14) 

There is a problem with this argument: we are using fixed-spin boundary 
conditions, not cyclic ones. However, provided the conditions are commen- 
surate with the ground state of the system, we expect (13. 14j) to be true in 
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the limit when the lattice is large and only a finite number of the spins 
in p, p' differ from the ground-state values. This is how we construct the 
infinite-dimensional matrices, so we believe (13.141) to be true in this limit. 
For i odd, it follows that the matrices 

AA i+l , BiB i+1 , A^Fi B i+1 , J B i G i A i+1 (3.15) 

are independent of horizontal rapidities such as q, q' . Similarly,for i even, 
they are independent of vertical rapidities such as p,p'. 

Each matrix will depend only on the rapidities p, q of that particular, 
corner, row or column. For instance Bi will depend only on the rapidities of 
the lower-left corner. For each matrix we take all the p's to be that same, and 
all the g's, and write the matrix as a function of p, q, e.g. 

Bi = B l (p,q) . 

For i odd, each of the expressions (I3.15P is a product of matrices corre- 
sponding to blocks of C that occupy the same rows of C. The horizontal 
rapidity q must therefore be the same for each matrix within the expression, 
but the vertical p rapidities may be different. Thus for instance the above 
argument gives 

Bi(p,q)B 2 (p' ,q) = independent of q 

for all p,p',q. Fixing p' and assuming the Bi matrices to be invertible, it 
follows that there exist invertible matrices R(p), S(q) such that 

B 1 (p,q)=R(p)S(q) . 

This is a factorization property. All the corner matrices A^, Bi similarly 
factor and we find that 

Km = iximr^xU) , 

Bi(p,q) = [Ki-ifo)]" 1 ^*^) , (3.16) 

G-( P ,q) = [Ymr^ipdxtm , 

where pi = p, cji = q if % is odd; and pi + i = q~i, g, + i = pi for all i. Note that 
the matrices Ai , Bi are independent of p and q, while T { , Q i depend only 
on the single rapidity pi. 
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The choice of the X i? Yi is not unique: each can be pre- multiplied by dif- 
ferent constant invertible matrices Cj, D { and the A4, £>;, Qi then adjusted 
accordingly so as to ensure that (I3.16P remains satisfied. However, this does 
not change the eigenvalues of 



JVC = A\A^A%A\ > N = A1A2A3A4 . 

These matrices define and Z% and play an imortant role in the calcula- 
tions: for instance, the order parameters are 



) = ^c^TraceA^ / £ Trace .M a (3.17) 



where uj = e 2m l N and j — 1, . . . , N — 1. 

Substituting these forms into the equations (I2.9al) - (12.9dj) . all the Yi 



matrices cancel out, so the effect is to replace the matrices Ai,Bi,Fi, Gi by 
Ai, Bi, J-"i, Qi, where Ai,Bi are rapidity-independent and Ti,Qi depend only 
on the single rapidity pi. 



3.1 Properties and symmetries 

There are various properties and symmetries that relate the corner transfer 
matrices. For the Ising, self-dual Potts and chiral Potts models [13] it is true 
that 

W pp (n) = 1 , W pp (n) = 5 nfi (3.18) 

and there exist rotation and reflection transformations p — > Rp, p — > Sp, 
such that 

W pq (n) = W q , Rp (n) , W pq {n) = W q>Rp (-n) , (3.19) 

W pq (n) = W Sq ,s P {n) , W pq (n) = W Sq ,s P {-n) . (3.20) 

If we write £/j(n) more explicitly as Ui(p, q\n), it follows that Ui(q, Rp\n) = 
U i+1 (p,q\n) and 

4?(p,p) = Bt(p,p) = 1 forz = land3 , (3.21) 

A^Rp) = A a i+1 (p,q) G a i{q,Rp) = G a i+1 {p,q) , (3.22) 

and, for instance, 

Bi(Sq,Sp) = B 6 ^(p,q) T . (3.23) 
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Replacing p, q by q, Rp is simply equivalent to rotating the lattice through 
90°, i.e. incrementing the suffixes % by one. 

These equations impose many restrictions of the matrices Xi,Yi, A4, Bi, 
Ti, Qi- We shall not explore all of them here, but note that they do imply 
that there exists a matrix function X (p) such that we can choose 

Xlip) = Xl{p) = X\p) , X a 2 (p) = Xl{p) = (Al)' 1 X a (p) , (3.24) 

A a A a 

Al(p,q) = [x a (p)Y 1 X a (q) , A a 2 (p,q) = [x a (q)\ ~ % X\Rp) , (3.25) 

A a 3 (p,q) = [x a (Rp)]' 1 X a (Rq) , Al(p,q) = [X^Rq)]^ X\R 2 p) , 
where 



X (R 2 p) = M X (p) 



(3.26) 



Corresponding equations hold for the matrices Yi,Bi, B i: J\fi, with no super- 
fixes a. 



4 Ising model 

We have seen how the star-triangle relation implies that row-to-row transfer 
matrices commute and corner transfer matrices factor. To proceed further 
we need an important property, namely that for many models we can choose 
the rapidity variables p, q so that 

W pq (n), W pq {n) = functions only of q—p and n . (4.27) 

We call this the rapidity difference property. It is intimately connected with 
the fact that for such models the Boltzmann weights live on algebraic curves 
of genus or 1, so can be naturally be parametrized in terms of Jacobi elliptic 
functions. One had become so used to this property that it was taken for 
granted. It was the discovery of the chiral Potts model[13], which does not 
have this property, that made its importance obvious. 

The Ising model is a two-state model (i.e. iV = 2) of the type discussed 
above. Taking Wp^O) = 14-^(0) = 1, its other Boltzmann weights are[8] 

W pq (l) = e~ 2J = k' scd(fsT - u) , (4.28) 
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Wpq(l) = e~ 2J = k'scd{u) , (4.29) 

where 

u = u q — Up , (4.30) 
and the function scd(u) is defined by 

scd(w) = sn(u/2)/[cn(«/2)dn(«/2)] , 

snu, cn«, dnw being the Jacobi ellptic functions of argument u and modulus 
k. 

Then, setting k' = (1 — fc 2 ) 1//2 , we have 

sinh 2J = — , sinh2J = (4.31) 
cn k fc'sn u 

and sinh2Jsinh2J = The conjugate modulus k' is small at low tem- 

peratures, increasing to one at criticality. Within this range the model has 
ferromagnetic order. The Boltzmann weights satisfy (I3.19P if we define 

u Rp = u p + K , (4.32) 

they are real and positive provided u = u q — u p is real and 

< u q - u p < K , (4.33) 

K, K' being the usual complete elliptic integrals. Note that this restriction 
remains satisfied if we replace p, q by Rq, p, when u q — u p becomes K —u q +u p . 

The relation (I4.30p manifests the difference property. It has far-ranging 
consequences, in fact it enables us to calculate the spontaneous magnetiza- 
tion, i.e. the order parameter (13.171) with j = l.[U chapter 13] 

To see this, note that we write the first equation (I3.25P as 

Aliuq-up) = [X{up) a y l X{u q f . (4.34) 
Replacing p, q by q,r and multiplying, this gives 

A\(Uq - Up)A a i {u r - U q ) = A^Ur - Up) 

i.e. 

A*(u)Al(v) = Aliu + v) , (4.35) 



13 



for arbitrary u, v. Interchanging u, v, it follows that 

Al(u)Al{v) = Al(v)Al(u) , (4.36) 

so A^iu), A\{v) commute, for all u, v. 

There is therefore a similarity transformation, independent of u, that 
diagonalises A 1 (u). Let us assume this transformation has been applied and 
take A^(u) to be diagonal. Then it follows from (14.351) that each diagonal 
element of A^(u) is an exponential function of u, i.e. 

= exp(-niw) , (4.37) 

where is some number. 

Allow u = tig — Up to be an arbitrary complex number. Since 

scd(u + 2xK') = - scd(w) , 

it follows from (14.281) that incrementing u by 2iK' merely negates W pq (l) 
and Wpq{l). From series expansions, we expect A x (u), and indeed all the 
functions we have defined in the thermodynamic limit of a large lattice, to be 
an analytic function of u within the extended physical domain < 9ft(tt) < K, 
which implies 

A a 1 (u + 2iK / ) = (-l^Aliu) . (4.38) 

It follows that 

m = 7^m ^ /(2K , ) , (4.39) 

where rrti is an even integer for a = 0, an odd integer for a = 1. 

If we take the low-temperature limit k' — > 0, keeping u in the vicinity of 
K/2, then to leading order in q' = Q~ 7rK / K ? 

Wpq (l) = e -™/^ , W pq (l) = e^-WPK') . (4.40) 

We find that A[ (u), as originally defined before (j2.4j) . is already diagonal 
to leading order. Let the spins to the left of and including a in Fig. [2] 
be ...,A2,Ai,a, where Ai is next to a, A2 is next to Ai, etc. Then the 
corresponding diagonal element of A^u) is indeed of the form (14.371) . (I4.39p . 
the index i being replaced by the spin set A, i.e. 

Aj(u) AlA = exp{-irum x /(2K')} , (4.41) 
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where 

m x = \ a - Ai| + 3|Ai - A 2 | + 5|A 2 - A 3 | H . (4.42) 

We have obtained this result in the low-temperature limit k' — > 1. However, 
for all k! the m\ must be integers and we do not expect them to change 
disco ntinuously as k' increases to unity. Hence we expect the formulae (14.41 [) , 
( I4.42p for the diagonal elements of A^(u) to be exact for < k' < 1. 

Set Mi = w 2 = u q — Up and u 2 = w 4 = K — u q + u p . Then from f!3.22p . 

A"(p,q) = Al{ Ui ) , (4.43) 

so these matrices, together with the product 

M a = A a 1 ( P ,q)A a 2 (p,q)A a 3 (p,q)Al(p,q) (4.44) 

are diagonalised by the same transformation that diagonalizes A^(u), and 

<a = , (4-45) 

where x — q' — e -71 "^/^'. 

These are also the elements of the diagonalized matrix A4 a \ we see that 
our calculation is consistent with the independence of on the rapidity 
variables p, q. 

It is now straightforward to calculate the spontaneous magnetization from 
(I3.17p . First set Ao = a and 

h = \\i-X l+ i\ , i>0 , (4.46) 

so li = if Aj = Aj+i, otherwise k = 1. Take Aj = for sufficiently large i, 
then to modulo 2, 

a = A + Ai + A 2 H . 

and 

tti\ = A + 3Ai + 5A 2 H . 

From flBTTTj) and (jQ5|l . 

((-l) a ) = ^(-l)V^ /J2 xmx ■ (4-47) 
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Here a, m\ are the linear expressions immediately above and the sums 
are over l Q , l\, l 2 , . . . = 0,1, without restriction. Each sum therefore factors 
into a product of individual sums over /g? ^2? • • ■■> giving 



((-in 



(l-x)(l-x 3 )(l-x 5 )--- 



(4.48) 



(1 +x)(l +x 3 )(l +x 5 ) ■ ■ ■ 



Using known elliptic function formulae .[T6l 8.197.4] and [U 15.1.4b] the RHS 



which is the famous result announced by Onsager in 1949 [23J and derived by 

Yang in 1952. [26] 

5 Summary 

We have expicitly shown how the star-triangle relation implies that the corner 
transfer matrices factor, causing the reduced matrix functions , Bi, J 7 ?, Ql 
to depend on at most one of the two rapidity variables. The rotation (and 
reflection) relations enable one to express the matrices, for % = 1, ... ,4, in 
terms of the single 2 = 1 case. 

Curiously, this by itself does not seem enough to enable one to calculate 
quantities such as the order parameters: if one also has the rapidity differ- 
ence property (14. 27ft , then the corner transfer matrices commute and can be 
simultaneously reduced to diagonal form. The eigenvalues are necessarily 
exponentials in u = u q — u p and are periodic of period 2iK'\ this means they 
are defined by certain integers m\, and these integers can be obtained from 
a low-temperature limit. The order parameters follow immediately. 

It still seems strange that the method should fail for a model such as 
the chiral Potts model, which does satisfy the star-triangle relation, so in 
that sense is "solvable". As yet the method has only yielded finite series 
expansions. [61 [7] 

The order parameters have now been obtained by a specialization of the 
Jimbo, Miwa, Nakayashiki method, [TS] - [UJ but it would be interesting to 
rescue the corner transfer matrix technique for this model. For the models 
with the difference property, diagonalizing the matrices leads directly to the 
Jacobi elliptic parametrization. For the chiral Potts model, is there any way 
of expressing the matrices that leads to an analogous parametrization? 



is &; 1 / 4 , so 



((-ir> 



£1/4 = (1 _ ^2)1/8 



(4.49) 
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